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Abstract 

We study the thermodynamic phases of a gas of spm-1/2 atoms in the Hartree-Fock ap- 
proximation. Our main result is that, for repulsive or weakly-attractive inter-component 
interaction strength, the superfluid and ferromagnetic phase transitions occur at the same 
temperature. For strongly-attractive inter-component interaction strength, however, the 
ferromagnetic phase transition occurs at a higher temperature than the superfluid phase 
transition. We also find that the presence of a condensate acts as an effective magnetic field 
that polarizes the normal cloud. We finally comment on the validity of the Hartree-Fock 
approximation in describing different phenomena in this system. 
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In recent years, studies of mult i- component Bose-Einstein Condensates (BEC) have re- 
vealed a variety of interesting phenomena that reflect qualitatively different ideas from spin- 
less condensates [HEllSllllElEllZllHlin!. Gases of spin-1/2 P, spin-1 |2l 13 El El El E] 

and spin-2 |H1 Ej atoms have been studied in considerable detail, both experimentally and 
theoretically. New physical phenomena include fragmented states [Sj, coreless vortices jB] 
and complex spin dynamics [UElEj- Truely spin-1/2 particles are fermions. However, a gas 
of bosonic atoms each of which can occupy two different internal states can be treated as a 
gas of spin-1/2 atoms that obey Bose statistics (note that in nonrelativistic systems, there 
is no a priori connection between spin and statistics). Examples of such gases are mixtures 
of ^^Rb atomic gases in two hyperfine states and spin-polarized Hydrogen E]] • 

It is well known that in a gas of noninteracting spin-1/2 bosonic atoms, the emergence of 
superfiuid order is accompanied by the emergence of ferromagnetic order jTU^j. A natural 
question that arises in that context is whether these two orders would remain connected if 
we include the effects of interatomic interactions. Perhaps a clearer way to pose the question 
is as follows: is it possible to have one of those two orders (i.e. superfiuid or ferromagnetic) 
without having the other? One can give a number of elementary arguments that hint one 
way or the other. On the one hand, in the noninteracting system the connection between 
the two orders follows from the symmetry of the wave function and is not related to any 
thermodynamic arguments (except for saying that the normal cloud is unpolarized). That 
would suggest a rather robust connection between superfiuid and ferromagnetic behaviours. 
On the other hand, if one considers the form of the order parameters or the symmetry of the 
Hamiltonian, there is no reason to believe that the two orders must be related. Keeping in 

mind certain caveats, the superfiuid order parameter can be chosen as {ip-^i^r)) + {tp^^r)) , 

whereas the ferromagnetic order parameter can be chosen as X]ij('^l(r)c?jjV'i(i"))) where aij 
is the vector of Pauli spin matrices, and the creation and annihilation operators used above 
will be explained when we give the Hamiltonian of the system below. The noninteracting 
Hamiltonian obeys SU(2) symmetry, whereas the interacting Hamiltonian obeys U(1)®U(1) 
symmetry. Furthermore, it has been shown recently by Yang and Li that the gound state of 
this system with no spin-dependent interactions is ferromagnetic, even if the interactions are 
strong enough to completely suppress superfluidity [T3|. The energies of low-lying excitations, 
however, also drops to zero for such strong interactions. It is therefore not obvious, at first 
sight, whether or not one order can exist without the other at finite temperature. A similar 
problem for a Bose gas of spin-1 atoms with ferromagnetic interactions has been considered 
by Gu and Klemm ^3]. One clear difference between spin-1/2 and spin-1 atoms is that 
with the latter one has a parameter that determines whether interatomic interactions are 
ferromagnetic or antiferromagnetic, even in a rot at ion- symmetric system. With spin-1/2 
atoms, apart from some obvious cases that we shall explain shortly, it is in general not 
immediately clear whether a certain set of interaction parameters describes ferromagnetic 
interatomic interactions. In fact, as explained above, a rotation-symmetric gas of spin-1/2 
atoms always exhibits ferromagnetic behaviour. 

In the present Paper, we shall try to obtain the answer to the above-posed question in 
the Hartree-Foclf (i.e. mean-field) approximation. We do that by calculating the free energy 
as a function of the macroscopic thermodynamic variables and minimizing the free energy 
with respect to those variables. Since similar approximations, as well as a number of sys- 
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tematic field-theoretic calculations, are not reliable in predicting the superfiuid transition 
temperature in a spinless system (see e.g. [ISIini, also see Appendix A), our results are nec- 
essarily plagued by the same type of unreliability, and they cannot be considered conclusive. 
The advantage of using the the Hartree-Fock approximation, however, is that whenever it 
gives correct results, it gives them with a simple physical explanation. At the end of our 
treatment, we shall give our intuitive assessment of which results we believe describe real 
physical phenomena and which results we believe are merely artifacts of the Hartree-Fock 
approximation. 

One rather trivial phenomenon in the context of spin-1/2 atoms occurs when the two 
internal states have different internal energies. In that case, the gas can be almost completely 
polarized at microkelvin temperatrues, assuming chemical equilibrium between the two spin 
states is reached. Although that phenomenon can be considered one form of ferromagnetic 
behaviour, we are not interested in it. Another ferromagnetic behaviour that we do not wish 
to consider here occurs when the two intra-component interaction strengths are sufficiently 
different. In that case it can, under certain conditions, be favourable for the majority of the 
atoms to occupy a single spin state. The system we shall consider in this Paper, which is also 
the most experimentally- relevant system of spin-1/2 atoms, is the one where the total spin 
along some axis is conserved, but the total spin perpendicular to that axis is not conserved. 
In the language of spin systems, one says that the longitudinal spin-relaxation time Ti is 
much longer than the time scale of performing the experiment, whereas the transverse spin- 
relaxation time T2 is shorter than the time scale of performing the experiment. We take the 
longitudinal axis to be the 2;-axis. Note that in the situation described above, and assuming 
that one starts with no net magnetization along the 2;-axis, the macroscopic magnetization 
expected for a condensate of noninteracting atoms lies entirely in the xy-plane. 

Let us take a Bose gas of spin-1/2 atoms in a three dimensional box. The Hamiltonian 
of the system can be expressed as: 

H = ifkin + Hiat, (1) 

where 

^ki„ = -^JdrY. ^Im'Mr), (2) 

o-=T,i 

^int = Jdr ^^i(r)V^:.(r)V^.'(r)^.(r), (3) 

o-,o-' = t,i 

m is the atomic mass, ^^^(r) is an operator that creates an atom at position r in spin state a 
(the spin states t and | correspond to the z-component of atomic spin being equal to ±1/2, 
respectively), V'o-(r) is its Hermitian conjugate, g^a' = 'inh'^aa-a' /m, and a^e^i is the scattering 
length between two atoms in spin states a and a'. We assume that the number of atoms 
in the a and a' spin states is individually conserved, whereas the total spin in the xy-plane 
is not conserved. Furthermore, we shall assume that there is an equal number of atoms in 
the two spin states t and | [ItJ- In order to avoid dealing with the possibility of mechanical 
collapse of the gas or phase separation of the two components, we shall assume that both 
fi-TT and are positive and that g'^^^ < g^^gu- 

We now derive an expression for the free energy of the above-described system as a func- 
tion of the macroscopic thermodynamic variables and use it to determine the state of the 
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system for a given temperature. The macroscopic variables are the number of atoms in the 
condensate No-, the total spin of the condensate S^, the number of atoms in the normal cloud 
Nn, and the total spin of the normal cloud Sat. We use the canonical ensemble, where the 
constraint of fixed total number of atoms No + N^ is imposed explicitly. Since we are dealing 
with a macroscopic system, we treat the above variables as classical variables. A condensate 
of spin-1/2 atoms is ferromagnetic as a result of Bose symmetry ^TJ^j. Its total spin Sq is 
equal to No/2 [i.e. S"^ = No/2[No/2 + l)]. On the other hand, the only constraint on the total 
spin of the normal cloud is that |SAr| < N]^/2. Therefore, we express the x, y and z compo- 
nents of the total spin of the condensate as {No/ 2 sin 9o cos ipo, No/2 sin 9o sin (fo, No/ 2 cos Oo) 
and those of the normal cloud as {NMSN±cosipN-iN]^SNL^va.^N, —No/2cos0o)- Note that 
with the above values of condensate and normal-cloud spins, we have taken into account our 
assumption that the 2;-component of the total spin must vanish. For any given set of values 
of the thermodynamic variables, the free energy is given by: 

F=(^kin) + (^int)-T5 (4) 

where T is the temperature and S is the entropy of the system. The entropy S and the 
expectation values are calculated by considering all the different microscopic configurations 
corresponding to the given values of the macroscopic variables. We first consider the first 
and third terms of the free energy, which we denote by Fideai- The condensate does not 
contribute to those terms in the free energy. Therefore we only need to evaluate -Fideai for 
a normal cloud of Nn atoms with total spin Sn- As we shall see below, the interaction 
energy is constant to leading order for all the different microscopic configurations in the 
Hartree-Fock thermodynamic ensemble for given Nj^ and S'at. Therefore, -Fideai is given by 
the same expression that it takes in the noninteracting system. The free energy of a spinless 
noninteracting uniform Bose gas is given by [TF 



-^ideal = (-f^kin) " TS (5) 

^5/2 
5^3/2 



N^ksTLz-'-^], (6) 



where /c^ is Boltzmann's constant and z is given by 

^3/2(2) = ^ivAi, (7) 
00 I 

9A^) ^Y.Tp (8) 



2-nh^ 
\ mkBT 



nj^ = Nj^/V, and V is the volume of the sample. -Fideai of the system at hand is given 
by the sum of two terms of the above form of -Fideai for two independent Bose gases, 
one with N^{l/2 + stv) atoms and the other with N]y{l/2 — s^) atoms, where sn = 



s%^ + N^cos^eo/4Nf,: 

-Fjdeal _ _j_ 

NisfkBT \2 "'V y""^ 5(3/2(2;+)/ ' V2 ^"J \ 5'3/2(^ 



ideal _^^,^__\(^^. g5/2M\ ^ fl „ \ f 95/2{Z-) \ ^ ^^^^ 
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where z± are given by 

9Z/2{Z±) =71^X1 (yJ^±SN^ . (11) 

Note that since z+ cannot exceed the value 1, must obey the condition n]^Xl{l/2 + s^) < 
5f3/2(l). Also note that Fideai is a monotonically increasing function of sjy- H we take the 
limit Sat —> 0, we find that: 

-a + 74 + ^4 + 0(4), (12) 



NNksT 
where 



1 ^5/2 (^) 

a = Inz . . , (13) 

^3/2(2) 

7 = (14) 
£/i/2(2;) 



2^|/2(^) / 3^2^/2(2^) ^-3/2(2;) 



^1/2 (^) ^V2(^) 



(15) 



and z is evaluated from Eq. ()11|) with sat = 0. When 2 = 0, 7 = 2 and = 4/3, and as ^ — > 1, 
7 oc const. Vl — z whereas 77 decreses slightly from the value 4/3 and remains finite. We now 
calculate the interaction energy in the Hartree-Fock approximation. In that approximation 
we assume that there is no coherence between states of different relative momentum of a pair 
of interacting atoms (i.e. (a^^^^ 0^.20.20^30-3 0^40-4) vanishes unless the momenta of the creation 
operators match those of the annihilation operators, not just the sum of the momenta). 
The two-particle correlation functions needed to evaluate the interaction energy can then be 
straightforwardly calculated to give: 

(V'}(r)^}(r)^t(r)V^T(r)) = n^cos^ ^ + li _ ^^^L cos ' + 2n,n^cos' % (l - ^cosft^) 

(V;i(r)V^}(r)V^;(r)V^l(r)) = sin^ ^ + ^ f 1 + — cos + 2n,n^sin2 ^ (l + ^cos^V) 

Z / \ fiN / Z \ riN J 

{'0j(r)?/'|(r)^/'|(r)^/'|(r)) = (-^-^ j + 44± + ^o^Af"SAf±sin6'oCos(v?o - V^at) (18) 

where n^, = No/V. By direct substitution of the above correlation functions, one can find 
the expression for the interaction energy (Hint) of the system. The state of the system can 
now be determined by minimizing the free energy: 

P g 

+ y , (19) 

where 
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fo 



-nl + gn\ + 2gnonN 



~( I \2 y+ 2 



(20) 



(21) 
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4 

5'TT + fl'ii + 9m 



(22) 
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5'TT ± 



(23) 
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Note that the fourth and fifth terms in Eq. ()19|) . among others, result from exhange- 
interactions. Those two terms describe scattering processes where an atom with spin \ 
and another with spin | exchange their momenta. In the expression for the free energy (Eq. 
IT!I|) . the second and third terms favour maximizing | cos^ol (i-e. having finite values of the 
z-component of the condensate and normal cloud), whereas the fifth and sixth terms favour 
taking cos 6o = 0. The second term can be understood quite intuitively as follows. Since the 
interaction term is magnified in the normal cloud as compared to the condensate (because 
of the exchange terms), the normal atoms tend to accumulate in the spin state with less- 
repulsive interactions. For example, if we take g-^^ > g^, we find that the normal cloud will 
have an excess of atoms in the | spin state, leaving the condensate with an excess of atoms in 
the I spin state. Since in this Paper we are not interested in that phenomenon, we eliminate 
it by taking g^^ = g^. The third term in Eq. ()19|) can also be understood by considering the 
absence of exchange interaction terms within the condensate. Since we are requiring that 
half the atoms have spin | and the other half have spin j, it is straightforward to see that 
in order to gain the greatest reduction in interaction energy, the condensate atoms tend to 
accumulate in the same spin state (either | or |). The competition between the third term 
and last two terms of Eq. (fT^ determines whether the condensate and normal cloud will 
have any net magnetization along the z-axis. Clearly, at low enough temperatures, the third 
term will win, and one will have a finite ^-component of the polarization of the condensate 
and normal cloud (note that the net polarization of the entire cloud must vanish, which is 
one of our main assumptions). That result is quite interesting in its own right. However, 
since the main question addressed in this Paper concerns the phase transitions, we focus our 
attention on those relatively-high temperatures. We minimize the free energy with respect 
to Uo, cos 6* and sn± at any given combination of the parameters g^i, g-^i and T to find the 
thermodynamic phases of the system. The results of a numerical calculation of the order of 
the phase transition are shown in Fig. 1. It is also worth making some analytical remarks 
about the bahaviour of the system. For clarity we address the following two cases separately 
(in the following remarks, we shall implicitly use a result that we found from our numerical 
calculation, namely that cos6 remains negligibly small close to the transition temperatures) : 
Case 1: g^^ > 0. If Uo ^ 0, the free energy is minimized by taking — '{)n = (i-e the 
polarizations of the condensate and normal cloud point in opposite directions) and a finite 
value of SAr_L. In the limit Uo/un -C 1, 



Note that above the BEC critical temperature Tc, i.e. when Hq = 0, the normal cloud is not 




(24) 
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polarized at all. Below Tc the physics can be understood in terms of the condensate acting 
as an effective magnetic field that partially polarizes a paramagnetic normal cloud. Since 
the superfiuid transition occurs at a temperature where 'jkBT ~ g^^n (see Appendix A), we 
find that: 

71 

nNSN± < Y (25) 

Therefore, the polarization of the normal cloud is smaller than that of the condensate, and 
the net polarization does not vanish. Note that if we calculate the exact expression for 
nNSN±, substitute it in Eq. (fT^ and minimize F with respect to Uo (keeping Uo + un fixed), 
we find that the global minimum of F jumps discontinuously from a point with rio = to 
a point with a finite value of Uo- That would suggest a first-order phase transition to the 
superfiuid phase, which is also ferromagnetic. However, the Hartree-Fock approximation 
predicts a first-order phase transition in a spinless Bose gas jTHj (also see Appendix A), and 
we therefore suspect that this result must be an artifact of the approximation. 
Case 2: < 0. In this case the free energy is minimized by taking ipo — = ^ (i-e the 
polarizations of the condensate and normal cloud point in the same direction). The term 
9u^'n^'n± favours a polarized normal cloud. That suggests that the gas might exhibit ferro- 
magnetic behaviour even if rio = 0, as we shall see shortly. Using the small sn± expression 
for the ideal-gas free energy (Eq. IT^ we find that the free energy is minimized by choosing 

= \9nK keTnl \ ^ ^^g) 

if 'jkBT — \g'\i\n]\r ^ l^f-i-Jno. As in case 1, the condensate acts as an effective magnetic field 
that polarizes the normal cloud. One can also immediately see that the normal gas exhibits 
ferromagnetic behaviour when jhsT — \gii\nN becomes negative. Note that since 7 = 
when T/T° < 1, the ferromagnetic phase transition must occur at a temperature higher 
than that of the ideal-gas BEG phase transition T°. Assuming that no = just below the 
ferromagnetic phase transition, we find that that transition occurs at a temperature where: 

2^^^^^kBTierro = \gn\^N, (27) 
9l/2[Z) 

where z is given by: 

?!m=(^]"\ (28) 



5'3/2(l) VTferro/ 

Just below the transition temperature, the spin grows as: 

= ^jiF^^En. (29) 

Note that if \g^i\ is smaller than a certain value, the ferromagnetic transition temperature 
can be smaller than the superfiuid transition temperature, which is also higher than T° for 
positive g-^i- In that case there would be a single (first-order) phase transition, just as in the 
case g^i > 0. However, as can be seen from Fig. 1 and as explained in Appendices A and 
B, there is a region in parameter space where the ferromagnetic transition temperature is 
higher than the superfiuid transition temperature. 

We now summarize our results and comment on them. We found that when gi^i > 0, there 
is a single transition to a phase that has both superfiuid and ferromagnetic order, and the 
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phase transition is first order. We believe that the first-order nature of the phase transition 
is an artifact of the Hartree-Fock approximation. However, since the physical mechanisms 
favouring an unpolarized gas (interaction energy) and those favouring an unpolarized normal 
gas (entropy) are both real physical mechanisms, we suspect that as soon as there is a 
finite fraction of atoms in the condensate, the induced polarization of the normal cloud 
will be smaller than the polarization of the condensate, and there will in fact be a single 
phase transition to a state with both superfluid and ferromagnetic order. In the case g^i < 
0, we found that it is possible, for large enough \g^i\i to have two phase transitions. At 
temperature Tferro, we found a second-order phase transition to a ferromagnetic phase with 
no superfluid order. We suspect that that result will persist even beyond the Hartree-Fock 
approximation, because that transition occurs above the superfluid transition temperature, 
and therefore fluctuations in the superfluid order parameter are expected to be negligible. 
At a lower termperature we found a first-order phase transition to a superfluid phase (with 
ferromagnetic order). As above, we believe that the superfluid phase transition will also 
be second-order. Based on similar elementary arguments alone, we cannot comment on the 
location of the boundary separating the single-transition and double-transition regimes. 

In conclusion, we have performed a Hartree-Fock (mean-field) calculation to study the 
phase transitions in a gas of spin-1/2 bosonic atoms. We found that it is possible to have fer- 
romagnetic order with no superfluid order, but not vice versa. We also found that the phase 
transition to the ferromagnetic non-superfluid phase is second-order, whereas any transition 
to a superfluid phase is first order. We suspect that the result of a nonsuperfluid ferromag- 
netic phase describes a real physical phenomenon. Judging from the Hartree-Fock results in 
the spinless case, however, we believe that in the real system, all phase transitions will turn 
out to be second order in nature. Due to the unreliability of the Hartree-Fock approximation 
near the superfluid transition temperature, further study is required to confirm or refute the 
results of this Paper. 

The author would like to thank T.-L. Ho for useful discussions and Q. Gu for useful 
correspondence about the results of Ref. ^1]. This work was supported by the National 
Science Foundation through Grant Nos. DMR-0071630 and DMR-0 109255 and by NASA 
through Grant Nos. NAG8-1441 and NAG8-1765. 
Appendix A: 

In this Appendix we show that the Hartree-Fock approximation predicts a first-order 
phase transition in a weakly-interacting spinless Bose gas ^2] (Here we are assuming repulsive 
interactions). Using similar arguments to those used in the main text, we find that the free 
energy density as a function of condensate number density Hq and normal cloud number 
density n^r is given by: 

^ = n^{\^z- - + g{n^ + n^)^ (30) 



where 



27rr 



3/2 



with the constraint that the right-hand side of Eq. (jM]) cannot be greater than (73/2(1)- The 
dependence of {F{no)—F(fi))/V on Uo (keeping the total density fixed) is shown schematically 
in Fig. Al. Clearly this calculation predicts a first-order phase transition at a temperature 
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Tc > T°, where T° is the transition temperature of the noninteracting system. By expanding 
F/V in powers of Uo, we find that, to leading order in na^, the shift in transition temperature 
is given by: 

= 1.08{naY'- (32) 



Note that as long as > 0, we find that z < 1 even below the transition temperature. Also 
note that when applying the above analysis to the situation discussed in this Paper, e.g. by 
taking = 0, some additional factors of 2 appear that lead to replacing the factor 1.08 in 
Eq. ((221) by 0.86, with that equation now relating (T^ - T°)/T° to {na^^^y/^. 
Appendix B: 

In this Appendix we calculate, to leading order the transition temperature to the ferro- 
magnetic nonsuperfluid phase. As explained in the main text, we treat only the case < 0, 
and we assume that \g^i\ is large enough that such a phase exists for a certain range of tem- 
peratures. We take Eq. (j2Z|), and use the asymptotic behaviour of the functions (73/2 (-z) and 
gi/2{z) as follows: 

hm 3V2{z) - 93/2{l) ^ ^gg^ 

z^l- y'l — Z 

which gives the asymptotic function: 

93/2iz) - (73/2(1) ^ -3.545v^r^. (34) 

The asymptotic limits of the functions gi/2{z), g^i/2{z), ... can be derived using the above 
approximation for (73/2 (-2) and the relation dgj{z)/dz = gj^i{z)/z. Using the asymptotic 
forms of (73/2(2;) and gi/2{z) in Eq. ^Tf^i. we find that 

^^^ = 1.84(n|an|f'' (35) 

c 

Comparing the results of Appendices A and B, one would expect the boundary between the 
two regions, i.e. those corresponding to a single and double phase transitions, to occur when 
I'^Til/'^TT ~ 0.47. The value of the ratio I'^til/'^TT ^"^^ actual boundary will be higher than 
0.47 if one takes into account the fact that a negative value of a-^i reduces Tc, in agreement 
with the results of our numerical calculation shown in Fig. 1. 

Figures 

Fig. 1: Regions of different behaviour in a gas of spin-1/2 atoms. The x and y axes are 
given by gi = g^^n/kBT° and (72 = 9u''^I^bT°. In region I we find a single phase transition 
to a superfluid ferromagnetic phase. In region II we find two phase transitions, with the 
ferromagnetic transition occurring at a higher temperature than the superfluid transition. 
Regions III and IV were not considered in this Paper because of possible instabilities: the 
former being susceptible to phase separation between the two spin species, and the latter 
susceptible to the whole cloud imploding because the net interatomic forces are attractive 
in that region. 

Fig. Al: Schematic plot of the free energy F{no) — F{0) of a spinless Bose gas as a 
function of condensate fraction rio for T > (solid line), T = (dashed line), > T > T° 
(dash-dotted line) T = T° (dash-double-dotted line). The dotted line is the x axis and serves 
as a guide to find the temperature of the first-order phase transition. As the temperature 
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is lowered, the point that minimizes F jumps discontinuously from to a finite value at 
temperature that is higher than the noninteracting transition temperature T°. 
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